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Our work in numerical modeling is focused on the use of the 
random vortex method to treat turbulent flow fields associated with 
combustion while flame fronts are considered as interfaces between 
reactants and products, propagating with the flow and at the same 
time advancing in the direction normal to themselves at a prescribed 
burning speed. The latter is associated with the generation of 
specific volume (the flame front acting, in effect, as the locus of 
volumetric sources) to account for the expansion of the flow field 
due to the exothermicity of the combustion process. The model was 
applied to the flow in a channel equipped with a rearward facing 
step. The results we obtained revealed the mechanism of the forma- 
tion of large scale turbulent structure in the wake of the step, 
while it showed the flame to stabilize on the outer edges of these 
eddies^ . 

In the course of this year we concentrated upon three topics: 

1. Study of the fundamental aspects of the modeling technique 

2. Application of the model to the formation of a turbulent jet 

3. Development of novel concepts bringing forth the aerodynamic 
properties of turbulent flames 

On the first subject, three efforts have been successfully pur- 
sued. First is the development of the random element method, pre- 
senting the application of the random walk technique to model dif- 
fusion of energy. The concept of vortex sheets and vortex blobs is 
extended for this purpose to heat transfer sheets or temperature jump 
elements, and blobs of internal energy, that carry thermal gradients 
from heat conducting walls to the interior of the field and redistribute 
them by diffusion. Table 1 presents the fundmental idea of random 
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walk modeling of diffusion, exploiting the analogy between the 
Green junction of the diffusion equation and the probability density 
junction of a Gaussian random variable. Figures 1 and 2 show a 
comparison between the exact solution and the numerical solution for 
the problem of heating an infinite solid by an isothermal wall at 
y = 0, in terms of the heat flux profiles and the temperature 
profiles. Figures 3 and 4 depict the temperature profiles for a 
finite solid bounded by two isothermal walls, and an isothermal and 
an adiabatic wall, respectively. The extention of the procedure to 
handle two-dimensional diffusion is presented in Table 2. A one- 
dimensional approximation close to the walls is employed to imple- 
ment isothermal boundary conditions, while energy elements are used 
in the interior to diffuse the energy there. Figures 5 and 6 display 
the temperature profiles along the diagonal of a corner and a square, 
respectively, as compared to the corresponding exact profiles. 

The second fundmental problem we have solved is concerned with the 
generation of vorticity by the interaction between the pressure field and 
density gradients — a mechanism which is of importance to both flames 
and buoyancy effects. The motion of the temperature jump elements, 
governed by. the algorithm of the random element method, generates ele- 
ments of vorticity in the interior of the field while additional 
elements are generated by the no-slip condition at the walls. Table 3 
describes how random walk is used in conjunction with the principle of 
time splitting to solve the system of equations that describes natural 
convection flow over a vertical infinite isothermal wall. Figure 7 

shows a comparison between our solution (thick lines) and a finite 

(3) 

difference solution for two values of the Prandtl number' . 

The third problem in this category is that of flame propagation. 

In our original model an interface advection and propagation algorithm, 
based on the assumption that the flame is a jump in density, was used to 
propagate the flame. In order to treat flame propagation as that of 
a reacting surface governed by a finite rate of reaction, the problem 
was recast in terms of a reaction-diffusion equation in temperature. 

The algorithm of the random element method was extended to solve this 
problem by adjusting the strength of the temperature jump elements as 
they move according to the rate of reaction as described in Table 4. 
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The algorithm was tested by solving problems by the use of finite 
difference and finite element methods, demonstrating its capability 
of calculating flame propagation with proper accuracy. Figures 8 and 
9 show a comparison between the exact solution, evaluated for 
f(T) = T(l-T), and the numerical computations when the integration 
of the reaction part of the equation is done using a first-order 
Euler scheme in the first case, and an exact integral in the second 
case, respectively^). 

On the second subject, the random vortex method was applied to 

the problem of the formation of planar, two-dimensional turbulent jet 

at high Reynolds numbers. Figure 10 describes schematically the 

elementary processes of the random vortex method and how they are 

implemented to solve the convection-diffusion equation. The results, 

expressed in terms of the development of the vorticity field, are 

presented in Fig. 11. They reveal the formation of large scale 

turbulent eddy structures on both sides of the jet with a potential 

core inside. Few jet widths downstream , the two layers start to 

interact and the flow becomes dominated by the pairing of eddies from 

both sides. The turbulent eddies grow by entraining the non-turbulent 

fluid, while their trajectories become more and more convoluted as a 

result of the interaction between positive and negative vorticity. 

These pictures display a remarkable resemblence to experimental photo- 
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graphs obtained by Dimotakis, et al . v ' across the plane of symmetry 
of an exi symmetric jet. 

The concept of flame aerodynamics has been developed on the basis 
of our numerical modeling studies, supported by experimental observa- 
tions of turbulent flames. The ultimate conclusion of these studies 
is the dominance of large scale eddies over the flow field, while the 
flame itself becomes established at the outer edges of these eddies 
and acts as a semi -permeable membrane encompassing the burnt gases. 

It is the interaction between. these eddies and the expansion associated 
with the exothermici ty of combustion process that produces the char- 
acteristic aerodynamic pattern that we are now studying^). 
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TABLE 1. FUNDAMENTAL IDEA 


Differential Equation 

3 2 4) 

3t 3y 2 

Boundary Conditions 

4>(y,0) = 6(0) ; <J>(+ °°,t) = 0 

Constraint 

rOO 

$ = / 4> dy = 1 

— 00 

Formal Solution 

cf) = — — exp[-(y/a)^] ; a = /2at 

/ 2iia 

Stochastic Solution 

Z 6$ i = 1 

y i (t + At) = y. (t) + n i 
E(n^ ) = 0 ; E(n.j 2 ) = 2a At 

Local Sampling 

♦ My-y,) 

Global Sampling 

$ = Z 6 H(y-y i ) 


Nomenclature 

H - Heavyside step function 
6 - Dirac delta function 

E - Expecte'd value 
n - Gaussian random variable 
a - Standard deviation 
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TABLE 2. A HYBRID SCHEME FOR TWO-DIMENSIONS 


DOMAIN 

BOUNDARY 

y < <5 $ 

INTERIOR 

y > « s 

Differential Equation 

3T _ 3 2 T 

9t a 3y 2 

— = aV 2 T 
3t 

Boundary Condition 

T = 1 

T = 6(y-6 s ) 

Diffusing Element 

6T i 

* e i 

Coupling 

6e i = ST.j * (y i -6 s ) 

«T f - ± 6e f /« s 

Stochastic Solution 

r. (t+At) = r^t) + n i 

(t+At) = r. (t) + n. 

Sampling 

Z 6e. 

T = 6 s -dx + E« T iH(y-y 1 -) 

T “ 6A, E 6e i 


$2 = 2a ; thickness of 1-D diffusion layer 
6A.j = area element 
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TABLE 3. PRESSURE-DENSITY INTERACTION 



n - number of temperature jump elements 
N = Number of vortex elements 
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TABLE 4. REACTION-DIFFUSION EQUATION 


DIFFERENTIAL EQUATION 
INITIAL CONDITION 
BOUNDARY CONDITIONS 

* J? { ) 

T(x,0) = T q (x) 

T(+°°,0) = 0 ; T(-°°,0) = I 

SOLUTION 

FRACTIONAL STEPS 

REACTION 

DIFFUSION 

T(nAt) = [R(At) D(At)] n T q 

ar= f < T > 

3T _ 3 2 T 

dt 3x 2 

DETERMINISTIC 

61. {t + At) = 6T.(t) + f(T.) At 

STOCHASTIC 

T = l 61. H(x- Xi ) 


x^ (t + At) = x i (t) + n i 


E[n] = 0 ; E[n 2 ] = 2 At 


256 






0 0.2 0 A 0.6 „ 0.8 1.0 

T 


Fig. 2 


257 









2.0 



Fig. 7 



Fig. 8 
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COMBUSTOR FLAME FLASHBACK 


James S. T'ien and Margaret Proctor 
Case Western Reserve University 
Cleveland, Ohio 44135 


Objective : Analytically model and conduct basic experimental tests to study 

the fundamentals of flame flashbacks in combustors. 

Background : Flashback has been a recurrent problem with the present Chrysler 

automotive gas turbine combustor. It can be a potential problem for 
other types and applications of advanced combustors as well. An improved 
understanding of the phenomenon would lead to improved design techniques 
to avoid its occurrence. 

Approach : It is proposed to model, test, and develop the fundamental 

conditions and process by which flashback occurs on combustor wall 
surfaces. An atmospheric rig consisting of a small dump combustor with a 
premixing channel would be required. Operation would be on a gaseous 
fuel. The primary variables would be inlet air and wall temperatures (up 
to 1500° F), boundary layer thickness, gas stream velocities, and 
controlled pressure disturbance level in the combustor. Both 
steady-state and transient tests will be conducted. 

Status : A steady-state flame flashback model is in existence. A transient 

model is being considered. The experimental rig is being designed. The 
experiment will be performed at NASA Lewis Research Center by a Case 
Western Reserve University student. 
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